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Abstract. 

We present iron abundance measurements, based on high resolution spectroscopy, and accu- 
rate distance determinations, based on near infrared photometry, for 34 Galactic Cepheids. 
The new data are used to constrain the Galactic iron abundance gradient in the outer disk, 
namely from 10 to 14 kpc. We confirm the flattening of the gradient toward the outer disk. 
In this region we also found an increase in the metallicity dispersion. Current data do not 
support the occurrence of a jump in the metallicity gradient for Galactocentric distances of 
the order of 10-12 kpc. 
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1. Introduction 

Galactic abundance gradients are a funda- 
mental input for chemodynamical evolutionary 
models, since they are the observables typi- 
cally adopted to validate predictions. Among 
the different tracers used to determine the 
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* Based on observations obtained with 
ESPADONS at the Canada-France-Hawaii 
Telescope (CFHT) and on observations collected 
with FEROS at the ESO/MPI 2.2 m. telescope at 
ESO, La Silla. 



Galactic gradients: HII regions. Open Clusters, 
O/B-type stars. Planetary Nebulae, we chose 
the Cepheids. The Cepheids present several ad- 
vantages when compared with other tracers: i) 
they are excellent distance indicators, which 
is a strong positive feature to provide accu- 
rate Galactocentric distances; ii) they are also 
bright enough to allow the study of the gradient 
over a large range of Galactocentric distances; 
Hi) they present a large set of well defined ab- 
sorption lines, therefore, accurate and precise 
abundances of many heavy elements can be 
provided. 
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Abundance gradients were discovered first 
in six ex t ernal galaxies from HII regions 
bv ISearld (Il971h but their occurrence in the 
Galaxy was a controversial issue. Indeed, some 
investig ations were in f a vor o f Galac t ic gra - 
dients jP'Odorico et all 119761; iJanesI 119791) . 
while others found no correlation between 
metallicity and distanc e (^legg & Bell 1973; 
iJennens & He"ife3ll975h . 

Even though the presence of Galactic 
abundance gradients seems widely accepted, 
the empirical estimates of the slopes are 
still lively debate d. By u sing HII regions, 
IVilchez & EstebanI (Il996h found a slope 
of -0.02 dex kpc ' , but different authors 
using the same tracers suggest slopes ra nging 
from -0.039 dex kpc ' (Deharv eng et al 



200C) to -0.065 dex kpc 
19971) . On the other hand 
slopes 



(iAfflerbach et al 



the use of B-type 
ranging from -0.042 
I2004 J to -0.07 
119981) . The 



stars gives 

dex kpc ' CDaflon et al 
dex kpc ' dGummersbach et al 
slopes based on Planetary Nebulae 
from -0.05 dex kpc 

the 



range 

(Costa et all l2004l) 
0.06 dex kpc ' (Mac iel et al.l 1 19991) to 
lack of a Galactic metallicity gradient 
dStanghelUni et al.l2006h . The slopes based on 
old Open Clusters still show a large spread. By 
adopting a sample of 40 clusters distributed 
between the solar circle and Rq - 14 kpc, 
iFriel et al.1 (l2002l) found a slo pe of -0.06 
dex kp c"' . More recently, ICarraro et al.l 
(1200 7') using new accurate metal abundances 
for five old open clusters located in the outer 
disk together with the sample adopted by 
Friel et al. (2002) found a much shallower 
global iron gradient, namely -0.018 + 0.021 
dex kpc ' . The global iron slopes based 
on Cepheids are very homogeneous, and 
indeed dating back to the first estimates by 
iHarrisI (1 19811 11984 . who found a slope of 
-0.07 dex kpc"' , the more recent estimates 
provide slopes ranging f rom -0.06 dex kpc"' 
(lAndrievskv et al."2002a''b'Q; 'Luck et al."2003t 
[Andrievsky et al. 2004; Lucketal. 2006) to 
-0.07 dex kpc"' (ILemasle et al.ll2007i) . 

Concerning the shape of the gradient, 
the situation is even more complicated. The 
hypothesis of a linear gradient is very dis- 
puted. Several authors favor a flattening of the 



gradie nt beyond 10-12 k pc: Vilchez & EstebanI 
dl996h (HII regions), 'Twar og et al.l (|l997h 
(Open Clus ters), lA ndrievsk y et al.1 d2004) 
(Cepheids), 'Costa etaU (12004 *) (Planetary 
Nebulae). This flatt ening is also well re pro- 



duced by models ( Cescutti et al.l 



2007h . In 



recent paper, lYong et aL (2006) brought 



another feature: the flattening may occur with 
two basement values, the first one at -0.5 
dex, a lower value than previous studies, 
and the second one at -0.8 dex, for which 
they suspected the possibility of a merger 
event. On the other hand, independent in- 
vestigations do not show a clear flattening 
toward the outer dis k, like RoUeston et al.l 
d2000l) (0,B stars) and lDeharveng et al.1 d20()0h 
(HII regions). Morever, a change in the 
slope in the direction of the inner disk was 
suggested by lAndrievskv et al.l (l2002bl) in a 
study based on a limited sample of Cepheids. 
In this context it is worth mentioning that 
the different tracers adopted to estimate the 
global iron gradient present a handful of 
objects in the outer disc, i.e. at Re > 12 
kpc. Finally, we mention that it has also been 
suggested by Twaro g et al. (1997), usi n g Ope n 
Clusters, and by lAndrievskv et al.l (|2004 . 
using Cepheids, a jump in the metallicity 
gradient at Rc ~ 10-12 kpc, with the iron 
abundance sharply decreasing by « -0.2 dex. 
Here we investigate the shape of the gradient 
in the outer disk and focus our attention on the 
possible jump in metallicity. 



2. Observations and Methods 

Our data sample includes high resolution 
spectra of 34 Galactic Cepheids. A large 
fraction of them (28) were collected with 
ESPADONS at CFHT, whereas six were 
collected wifli FEROS at 2 2m ESO/MPG 
telescope dKaufer et al.l Il999l) . Spectra were 
reduced using either the FEROS package 
within MIDAS or the Libre-ESpRIT software 
dDonati et al. 1997, 2006). The data analysis 
was already described in lLemasle et al.l (120071) . 
Suffice it here to briefly mention the main steps 
of the reduction strategy. The accurate deter- 
mination of effective temperatures is a criti- 
cal point in the abundance determination. They 



Lemasle et al.: Cepheids as tracers of the Galactic disk metalhcity gradient. 



3 



are only spectroscopically determined, by us- 
ing the me thod of line depth rati o s (LD R), de- 
scribed in iKovtvukh & Gorloval (l2000h . As it 
relies only on spectroscopy, it is independent 
of interstellar reddening. These authors pro- 
posed 32 analytical relations for determining 
Tetf from LDR of weak, neutral metallic lines. 

The other atmospheric parameters (surface 
gravity log g , micro turbulent velocity v,) are 
determined by imposing the ionization balance 
between Fel and Fell with the help of curves 
of growth: iterations on log g and v, are 
repeated until the best match with the same 
curve of growth is reached. Abundances are 
then calculated with an atmo spheric model 
from lEdvardsson et al.l (Il993h based on the 
following assumptions: parallel plane stratifi- 
cation, hydrostatic equilibrium and LTE. We 
adopted the s caled-solar chemical a bundances 
suggested by iGrevesse et al ] (ll996h . The final 
internal accuracy in the abundance determina- 
tion is of the order of 0. 1 dex. 

Absolute distances were estimated using 
near infrared photometry {J, H, K-bwA) from 
the 2MASS catalog. The mean magnitude of 
these objects was estimated using the tem- 
plate l ight curves provided by ISoszvnski et al.l 
(l2005h . together with the V-band amplitude 
and the epoch of maximum available in the 
liter ature. For ST Tau a nd AD Gem the data 
from lBarnes et al.l d 19971) were also used, while 
for H W Pup the data from ISchechter et al.l 
( 1199 2') have been joined to the 2MASS ones. 
As templates are only for fundamental modes 
pulsators, the mean near infrared magnitudes 
for first overtone pulsators are only the sin- 
gle 2MASS measurement. The Cepheid ab- 
solute magnitudes were computed using the 
near infrared peri od-luminosity relatio ns re- 
cently provided by iPersson et al.l (l2004l) . with 
an LMC distance modulus of 18.50 mag. 
Pulsation periods and reddenings are from the 
Fernie's database (iFernie et al.lll995h and pe- 
riod of first overtone pulsators were funda- 
mentalized. We assumed a Galactocentric dis- 
tance of 8.5 kpc for the sun (iFeast & Whitelocfl 
Il997h . 



3. Results and final remarks 

Current iron abundances when compared with 
previous measurements, mainly from the sys- 
tematic investigation by Andrievsky et al., 
show a very good agreement. For 17 stars 
the abundance difference is very small (<0.1 
dex) and for 6 stars it is ~ 0.1 dex. We 
cannot confirm previous determinations only 
for 4 stars, for which discrepancies can reach 
0.3 dex. To improve the sampling along the 
Galactic disk, we added to our new data 
30 Cepheids from iLemasle et all (l2007h . 52 
Cepheids from Andri evsky ' s sample and 11 
Cepheids from M ottini et al.1 (12007,) . for which 
both metallicities and infrared photometry 
were available. The complete sample includes 
127 Cepheids with distances based on homo- 
geneous infrared photometry. Among them 27 
have Galactocentric distances in the 10-12 kpc 
range and 13 are located beyond 12 kpc. 

We first estimated the Galactic gradient for 
the whole Cepheid sample. Data plotted in Fig. 
[Ubring forward the following findings: 

• The iron gradient shows no evidence of a 
gap over the Galactocentric distances covered 
by the current sample. 

• The iron gradient is flattening in the outer 
disk. 

• The spread in metallicity is larger in the outer 
disk, possibly due to local inhomogeneities. 

The slope of the linear iron gradient 
over the entire sample is shallower, -0.058 
dex kpc~' , th an pr eviously estimated by 
iLemasle et al.l (|2007|) . but the new sam- 
ple extends toward the outer disk, where 
the flattening occurs. This slope is in very 
good agreement with p revious determinations 
from Ei el et al. (l2002h (Open Clusters) and 
iLuck e7al.l (120061) (Cepheids). The flattening 
of the gradient in the outer disk becomes even 
more evident in Fig.|2] We estimated the slope 
by using Cepheids located between 10 and 15 
kpc and we found -0.019 dex kpc ' . This 
estimate is in excellent agreemen t with the 
slope found by ICarraro et al ] (l2007h using old 
Open Clusters, i.e. -0.018 + 0.02 dex kpc ' . 
Moreover and even more importantly, the two 
independent estimates suggest, within the er- 
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rors, a very similar metallicity value in the 
outer disc, namely {FelH} - -0.3, -0.4. 

Several reasons can be invoked to explain 
the lack of the metallicity gap i n our sam- 
ple. Th e first one was suggested by 'Luck et alJ 
(l2006h . who pointed out that their sample, as 
well as the one by Twaro g et al.l ( Il997h was 
gathering stars lying approximately at the same 
Galactic longitude, respectively 190-250 and 
130-260 degrees. Moreover, the uncertainties 
on the distance determination increase with the 
star distances (especially when the distance is 
estimated using optical photometry) and could 
introduce artifact in the distribution. The use 
of IR photometry and robust primary distance 
indicators presents several undisputed advan- 
tages. Note that the outer disk in strongly un- 
dersampled, therefore, the selection of the tar- 
gets among the few available Cepheids can in- 
troduce a bias that can only be removed by in- 
creasing the sample size in this region. Finally, 
we note that the lack of a gap in our iron gra- 
dient estimate might be due to the restricted 
range of longitudes covered by our sample. 
This relevant point deserves further investiga- 
tions. 



4. Summary 

High resolution spectra and accurate distance 
determinations for 34 Cepheids provided the 
opportunity to investigate the iron gradient 
across the Galactic disc. Preliminary results in- 
dicate that: 

• The Galactic iron gradient presents a linear 
trend and the slope is k -0.06 dex kpc ' . 

• The shape of the gradient is more accurately 
described by a bimodal distribution, with a 
higher slope in the central region and a flatten- 
ing of the gradient in the outer disk. In this re- 
gion, the spread in metallicity is higher, possi- 
bly due to local inhomogeneities. A third zone 
can be considered in the inner disk, with a 
break in the slope near 7 kpc, but this region 
was not explored in our study. 

• Current data show no evidence of a jump in 
the iron gradient for Rq ~ \0 — \2 kpc. 



SLOPE: -0.058 
I -sic^o; 0,004 



Fig. 1. Galactic radial abundance gradient. The 
solid line shows the linear regression whose 
slope and \ - cr error are labeled. Diamonds 
mark our data, while asterisks show data from 
Andrievsky et al. and plus signs display data 
from Mottini et al. 
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Fig. 2. Same as Fig.l, but for Cepheids with 
Galactocentric distances ranging from 10 to 15 
kpc range. 
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